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ABSTRACT 


This paper describes a proposed pulse-frequency mod- 
ulation digital communication system for satellite applica- 
tions. The heart of the proposed system is a maximum- 
likelihood detector of the spectrum-analysis type discussed 
by Darlington. It is shown that the proposed technique is in 
principle an ideal M-ary system. A general description is 
followed by an analysis of anticipated departures from the 
theoretic system— in particular, waveform distortion due to 
finite bandwidth, slope mismatch, and phase and amplitude 
ripple. General results are derived, from which permissible 
tolerances for departures from the theoretical system may 
be computed. 
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WAVEFORM DISTORTION IN A MAXIMUM-LIKELIHOOD 
DARLINGTON DETECTOR FOR M-ARY PFM SIGNALS 


by 

Paul J. Heffernan 

Goddard Space Flight Center 


INTRODUCTION 

This paper presents an analysis made at GSFC in 1966 to support a feasibility study for a 
time-division-multiplex system of multiple -access communications between pairs of small ground 
terminals via a subsynchronous earth-orbiting satellite. 

The modulation-demodulation equipment has not been instrumented beyond breadboarding 
several critical elements; timing and synchronization hardware has been developed to engineering- 
model level and tested in satellite links. 

While a brief discussion of communication system operation and anticipated performance is 
included for purposes of completeness, the main purpose of this paper is to predict the effects of 
amplitude and phase anomalies in a particular time -invariant linear system using classical Fourier 
integral techniques. 


OVERALL SYSTEM OPERATION 
General 

This communication system furnishes duplex 4 kHz voice or data channels simultaneously to 
several pairs of users via an earth-orbiting satellite relay. All communication is in real time; 
multiple access is accomplished by time -division- multiplex (TDM). Multiple-access satellite com- 
munications have been treated extensively in the literature (References 1, 2, and 3), so it is enough 
to state here that of the three basic signal- multiplexing techniques available (time -division, 
frequency-division, and common-spectrum), TDM in theory provides the greatest channel capacity 
per user in situations where satellite transmitter power is the parameter limiting performance 
(References 4 and 5). 

The modulation technique selected is 127-level quantized pulse-frequency modulation (PFM). 
Signal detection is by "spectrum analysis" pulse-frequency modulation to pulse -position modulation 
(PFM-PPM) conversion. This scheme was first discussed by Darlington in connection with threshold 
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extension for wideband FM signals (Reference 6). The 127-PFM signaling waveforms are pair- 
wise orthogonal, and the Darlington demodulation process is essentially non-coherent matched- 
filter maximum-likelihood detection. Hence, the signaling scheme is in principle an ideal M-ary 
system (Reference 7). That is to say, the operations by which the receiver guesses which of the M 
signals was sent afford the least probability of error possible with non-coherent reception (in the 
absence of subsequent error detection and/or correction). 


Timing and Synchronization 

This discussion of timing and synchronization is included in order to define the communication 
system under consideration and is not a description in depth. The interested reader is referred to 
Reference 8 for a complete description of system timing and synchronization subsystems and to 
Reference 9 for the results of timing and synchronization subsystem performance in closed-loop 
satellite tests. 


Up to ten ground stations may participate in the system; only one station at a time uses the 
satellite. Hence, the system provides a maximum of ten simplex links or (when pairs of ground 
stations communicate with each other) up to five duplex links. 


The system centers around a timing frame that contains one master synchronization burst and 
ten information bursts, one for each simplex link. There are 8000 frames per second, so that no 
storage or buffering is required for transmission of voice or data signals bandlimited to 4 kHz 
(assuming periodic Nyquist sampling) as long as only one such 4 kHz signal is to be transmitted 
per simplex link. (See Figure 1.) 

Of the participating ground stations, one is called the master station; the rest are called slave 
stations. The master station controls all system timing. In each frame, the master station trans- 
mits the master synchronization burst from which all slave stations derive their timing. This 
burst is pre-programmed to be an rf pulse modulated in frequency at 800.00 kHz and repeated at 
intervals of 125.00 microseconds where these values hold at the satellite. Because the satellite is 
in motion with respect to master and slave stations, Doppler effects must be considered in signal 
pre-programming. 
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Figure 1 —Format diagram. 



Each slave station monitors the master synchronization burst upon its transmission back from 
the satellite to derive slave local timing. Because of the Doppler effect, local timing will differ 
somewhat (in any case less than a small fraction of 1 percent) from the exact 125.00 microsecond 
frame rate at the satellite; in general, local timing will differ slightly for each ground station in 
the system. 

Once local framing is derived, intraframe timing must be established to prevent ground- station 
transmissions from overlapping upon arrival at the satellite. Intraframe timing is required also 
for operation of the signal demodulation equipment at each ground station— see below. 

As already noted, techniques and hardware for establishing and maintaining these timing and 
synchronization functions have been developed, proven, and documented; further details of TDM 
system timing and synchronization are not given here. The proper subject of discussion is a point- 
to-point communication system in which 4 kHz signals are transmitted as sampled data; sampling 
intervals and instants being known with arbitrary precision and accuracy at both ends of the link, 
while rf carrier phase is unknown. 


Pulse-frequency Modulation and Darlington Detection 

PFM has been used with considerable success in numerous telemetry applications (Refer- 
ence 10), but few attempts have been made to use this efficient M-ary signaling technique in com- 
munications systems for voice and data transmission. 


The information signal to be transmitted is assumed to be an analog waveform bandlimited to 
4 kHz and sampled periodically at the Nyquist rate; the clamped samples drive a quantized voltage- 
controlled oscillator to produce a train of rather low duty-cycle FM pulses as shown in 
Figure 2. 


PFM 

BURST 



■125 h SEC 


125 jjSEC 


Figure 2— PFM signals. 
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A single pulse is of the form 


S N (t) = cos 2rr (f c +Sf)t , -T 0 <t <T 0 


( 1 ) 


where 

f c = reference carrier frequency 

Sf = 100 kHz (basic increment of discrete frequency modulation) 
N = 0, ±1, ±2, • • *, ±64 
T 0 = 5 /^sec 

The interval between pulses is 


T s - 125 m sec 

and corresponds to the Nyquist rate for 4 kHz signals. The discrete frequency modulation NSf — 
constant during any one pulse but changing from\mlse to pulse in accordance with the amplitude 
samples of the information signal— is selected to be an integral multiple of the reciprocal of pulse 
width 2T 0 . As a consequence and despite their overlapping spectra, the 127 signaling waveforms 
are pairwise orthogonal over the signaling interval (Reference 11). 

The amplitude resolution of the system is identical with that of seven-bit pulse-code modula- 
tion and, as with PCM, logarithmic or piecewise linear encoding would normally be used to reduce 
the effects of quantization distortion (Reference 12). 

The optimum non-coherent detector for these PFM waveforms may be visualized as a bank of 
127 matched filters whose outputs are envelope-detected and sampled at the anticipated instant of 
match to select the filter with the greatest response (Reference 13). The filter matched to the 
Nth signaling waveform s N ( t ) would have a system function h n (f) given by the complex conjugate 
of the spectrum of s N ( t ) : 


*C(f) 


f S N ( t ) c~ i27Tft dt - 

•'-a 


( 2 ) 


Equivalently, the impulse response h N (t) of the Nth matched filter would be the S N (t) waveform 
run backward in time: 


h N (t) = S N (-t) = S„(t) . 


(3) 
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Figure 3 is a block diagram of this PFM demodulator. Such systems have been instrumented 
for small alphabets (Reference 14); but for large alphabets (e.g., M = 127) the hardware would be 
unwieldy and expensive. 

The alternate detection scheme shown in Figure 4 is a slight variant of a technique proposed 
by Darlington and referenced above. The 127 matched filters are replaced by a single circuit con- 
sisting of a function generator, a mixer, a dispersive delay line, and a tapped non-dispersive delay 
line. The operation of the detector may be described as follows. 

The signal to be detected is S N (t) as given by Equation 1. It is mixed with f(t), the out- 
put of the function generator. This function generator waveform is a linear FM pulse de- 
scribed by 


f(t) 


COS 277 


t + 



- T < t < T 

1 0 _ L _ 1 Q 


(4) 


where f r is an unspecified reference frequency. The instantaneous frequency f(t) is given by 
the time derivative of the phase and is 

f i (O ■= f r + Mt , -T 0 <t<T 0 (5) 

where ^ = 1.28 x 10 12 Hz/second. 

The instantaneous frequency is thus a linear function of time (see Figure 5). After mixing and 
removal of out-of-band spurious responses, the signal driving the dispersive delay line is of the 
form 


x n (O 


COS 277 


(f 0 ’NSf)t 



- T < t < T 
1 o - L - A 0 


(6) 


where f 0 - f c + f r is the carrier frequency of the pulse, and both the discrete frequency modu- 
lation NS f and the linear frequency modulation fut 2 /2 are present. 

The response of the dispersive delay line to the excitation x N (t) can be determined using the 
Fourier integral. The complex system function of the delay line is modeled (for present purposes) 
as 


H(f ) = A(f ) e^ (f > , ( 7 ) 

where the amplitude characteristic A(f) is that of an all-pass structure, 
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Figure 4 — Darlington PFM detector. 


>“ 

o 

z 

LU 

ZD 

o 

uu 

o' 



t, TIME 


Figure 3— Filter bank PFM detector. 


Figure 5— Instantaneous frequency of f(t) vs time. 


A(f) r 1 , -oo< f <oc 


(8) 


and the phase characteristic is 


tf(f) = J ( f " f o) 2 - ^ 


( 9 ) 


The time delay vs frequency for this structure is 


T(f ) 


6 


- f 




+ 


( 10 ) 










Parameter 1 0 represents an arbitrary 
constant delay; the f 0 - f/^ term represents a 
delay that decreases linearly with increasing 
(positive) frequency. (See Figure 6.) Struc- 
tures or filters of this type are usually termed 
dispersive, because they delay the spectral 
components of an arbitrary excitation non- 
uniformly and hence tend to spread out or dis- 
perse the signal in time. 

This need not always be true. For in- 
stance, if the input is as given in Equation 6, 
the unit amplitude rectangular -envelope input 
pulse is compressed into a narrow, peaked 
(sin x )/ x pulse, where the time to peak is de- 
termined by the discrete frequency -modulation NSf 



Figure 6— Time delay vs frequency 
of dispersive delay line. 

and the overall network delay t 0 . This 


is evident from an expression for network response. Let 


y N (t) = 3^(0 ' x(t ) * h( t ) 


(id 


where 


h(t) 


r 


A(f)e 


i<£(f) _i27Tft 


df 


( 12 ) 


Taking A(f ) equal to unity, corresponding to the postulated all-pass amplitude characteristic, 
and with <£(f) given by Equation 9, the response is a rather complicated FM waveform whose en- 
velope is given by 


y N (t) 




sin 277 T q [ M (t - t Q ) + NS fj 
[/*(t- t 0 )+NSf]' _ 


( 13 ) 


This result is derived in Appendix A. Before studying the properties of this response, we 
simplify Equation 13 as follows: The overall filter delay t 0 is significant because it determines 
how far structures may be approximated at the frequency of interest. But— for purposes of signal 
analysis— it is an arbitrary parameter that can be set equal to zero. The response is then 


y N (t) 




sin 277 Tq [/j.t + NSf ] 
2 tt T 0 [fj. t + NSf] 


(14) 
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Input and simplified output waveforms are 
shown in Figure 7 for the case N = 0. Equa- 
tions 6 and 14 and Figure 7 indicate that an ef- 
fective pulse compression has taken place and 
that the discrete PFM has been converted to a 
discrete PPM. 

Associated with and characterizing the 
pulse compression phenomenon is a signal proc- 
essing gain P which may be defined as the ratio 
of output to input pulse peak power: 

r = 4u T 0 2 . (15) 

In the present system where 

IX = 1.28 X 10 12 Hz 

and 

T Q = 5 fisec , 

we have 

r = 128 , i.e. , 21.0 db . (16) 

This signal-processing gain* is a fundamental figure of merit for the general pulse-compression 
phenomenon because r,as defined above, is a measure of the pulse-width compression ratio. With 
reference to Figure 7, output-pulse width is defined in terms of its (approximate) half-power points; 
the ratio of input to output pulse widths is 

2T o + 2£T 0 = 4 ^ T o 2 = r * 

r may also be considered a time -bandwidth product, that is, the product of the duration and 
total frequency sweep of the input rectangular-envelope frequency-modulated pulse: 

(2T 0 )(2mT 0 ) = VT 0 a = r. 

♦Also termed dispersion factor in the literature of pulse-compression radar. 



-5-4-3-2-10 I 2 3 4 5 

2fjT 0 t 


Figure 7— Input and output waveforms for 
the Darlington PFM detector. 
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The PFM-to-PPM conversion process is the heart of the Darlington demodulator. For N / 0 
the y N ( t ) waveform peaks at 


t 


NS f 


N = ±1, ±2, * • • , ±64 


and for any N, the nulls of y N (t) occur at 


(17) 


NSf M 
1 M 1 2 ^T 0 


M = ±1, ±2, 


The system parameters have been selected so that 


(18) 


l 

S f - 9 'r 

Z1 0 

Therefore, the output waveforms y N (t) interlace peaks and nulls as shown in Figure 8. This quasi- 
orthogonality of the y N (t) is, of course, due to the orthogonality of the original signaling wave- 
forms and is the basis for the subsequent maximum -likelihood detection. 

The output of the dispersive delay line 
drives a non-dispersive tapped delay line. 

Taps are located at intervals of ( 2 ^T 0 )" 1 or 
about 78 nanoseconds corresponding to the 
pulse intervals indicated in Figure 8. Using 
techniques of the general type outlined in Ref- 
erence 14, the tap with the greatest response 
is selected to complete signal detection. 

Signal-to-noise Ratios and Error Rates 

It is stated above that the PFM communication technique using Darlington detection described 
here is essentially an ideal M-ary system of the type discussed by Turin. The theoretical perform- 
ance of such systems has been studied for some time and is well documented. References 7, 10, 

14, 15, 16, and 17 are especially relevant. 

This paper is designed not to reproduce or summarize such results but to show that the present 
scheme is indeed an ideal M-ary system. To do this it is enough to show that the signal-to-noise 
ratio* at the output of the dispersive delay line in Figure 4 is just what would be obtained at the 
output of the matched filter containing the signal in Figure 3. 


Yo(t) 



2 hV 

Figure 8 — Quasi-orthogonality of the output waveforms. 


*In the usual sense of peak signal to mean-square noise power ratio under conditions of signal reception in wide-sense stationary white 
gaussian noise. 
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If the received signal is immersed in noise of (two-sided) spectral density r/2, the signal-to- 
noise ratio at the output of the Nth matched filter of Figure 3 at the instant of match is the "energy 
contrast ratio,” i.e., the ratio of signal energy to noise spectral density: 


2 

v 



s N 2 (t)dt 



(19) 


To compute the signal-to-noise ratio that can be obtained with the present system, we must 
identify the equivalent noise bandwidth of the overall demodulator. The derivation of Equation 13 
in Appendix A assumes an all-pass amplitude characteristic for the dispersive delay line and as- 
sociated transmission system, and leads to an output signal-to-noise ratio of zero; i.e., an infinite 
system-noise bandwidth. Actually it can be asserted that system bandwidths are finite, if only to 
ensure the realization of the dispersive delay line (or a suitable approximating structure) at the 
frequencies of interest. What is the noise bandwidth of the system? 

It is shown in Appendix B that the appropriate noise bandwidth of the dispersive delay line 
proper is twice the total frequency sweep 2 ju T 0 of any one of the x N (t). Provided that all prior 
structures are broadband, the system noise bandwidth is 4^T 0 , i.e., 25.6 MHz. For this case, we 
can compute the signal-to-noise ratio as follows. 

From Equation 13, the peak signal power is 


[y N 2 (t )] - 4 M t 0 2 

J max 

and from the above discussion the mean- square noise power is 


4 ? m t 0 . 


This leads to a signal-to-noise ratio of 


_4mT 0 2 _ Tq 

VT 0 M 


( 20 ) 


which is 3 db poorer than that attained in the matched filter case.* If we could narrow the system 
noise bandwidth to 2/; T 0 (as in chirp-radar systems for detection of a single known signal) the dis- 
persive delay-line detector would have the same signal-to-noise ratio as the matched filter, as 
argued in Reference 19. 


* Actually, the situation is somewhat more involved than the above argument indicates. The function generator waveform f(t) effectively 
gates the input noise process on and off and creates a non- stationary noise process at the output of the Darlington detector. Applying the 
analysis given in Reference 18 to the present system shows that the variance of the non-stationary output noise process has the con- 
stant value 4/xn T Q (for the present case where we assume no signal prefiltering) throughout the signal observation interval. 
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In fact, we can narrow the system noise bandwidth to a value closely approaching 2/i T 0 (al- 
though the dispersive delay line needs 4/xT 0 for effective PFM-to-PPM conversion of rectangular- 
envelope chirp signals) —prefilter the input signals S N (t) by suppressing all energy outside the 
band (considering now positive frequencies only) f c ± 6.4 MHz. This operation eliminates all 
out-of-band sidelobes of the sin x/x spectra of the S N (t) but need not affect in-band spectra in any 
way. As a consequence, there will be a new set of signals S N ' (t) which within the band have the 
same spectra as the original S N (t) but outside the band have different (i.e., highly attenuated) 
spectra. After mixing with f(t), the waveforms into the dispersive delay line will be of the form 
(using complex notation) 


*n (*) 


(t) 


i277 [( f o +NSf ) t+Jut 2 / 2 ] 


- T < t < T 

j. n — x. — x n 


(21) 


where e N ' (t) is a non -rectangular envelope function corresponding to the inverse Fourier trans- 
form of the spectrum of s N ' (t). If we ignore the explicit form of x N ' (t), the quadratic phase dis- 
persive delay line produces an output signal y N ' (t) whose envelope is the Fourier transform of 
e N ' (t) (see Appendix A). Within the signal observation interval, the y N ' (t) appears as pure sin x/x 
waveforms whose peaks and nulls interlace as discussed above; outside the observation interval, 
the response of the delay line is highly attenuated. 


This effect may be visualized by picturing the time response of the dispersive delay line as a 
scaled replica of the spectra of the prefiltered input waveforms; if prefiltering does not alter the 
spectra within the band f e ± 6.4 MHz, the time response of the dispersive delay line will not be 
altered or distorted within the observation interval; signal detection will be affected only in im- 
proving the signal-to -noise ratio of Equation 20 by an amount approaching 3 db. Hence, the present 
system, in the limit for perfect signal prefiltering and with a perfect Fourier-transforming dis- 
persive delay line, yields the same signal- to-noise ratio as a bank of perfect correlators or matched 
filters. 


WAVEFORM DISTORTION IN THE DARLINGTON DEMODULATOR 

Souces of Waveform Distortion 
and Philosophy of Analysis 

The derivation of the dispersive delay line response to the excitation x N (t) as given in Ap- 
pendix A assumes an infinite bandwidth system with unity all-pass amplitude characteristic and 
pure quadratic phase characteristic for all frequencies. A practicable system can, of course, only 
approximate these characteristics in some finite band, and it is to be expected that the actual out- 
put waveforms obtained will differ somewhat from those shown in Figures 7 and 8. 

Such departures from pure (sin x)/x waveforms are termed distortion and will, in general, 
produce two effects that may be considered separately. First, there will probably be some loss in 
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signal-processing gain r relative to the theoretical value of 128 or 21 db and second, the delicate 
quasi-orthogonality of the output waveforms will probably be destroyed to some extent, as the peaks 
and nulls may not interlace perfectly according to the scheme shown in Figure 8. 

This paper does not consider how far these distortions can be tolerated. It seeks general 
formulations characterizing losses in r and destruction of sidelobe interlacing in terms of depar- 
tures from ideal signals and networks considered in Appendix A. Emphasis is placed not on actual 
details of output signal waveshapes, but on loss in peak response (reflected in smaller output signal- 
to-noise ratios) and the rise, if any, of the (sin x)/x sidelobe level at the expected nulls, particularly 
those adjacent to the main response (see Figure 8). 

Potential sources of waveform distortion are finite system bandwidth, slope mismatch (Fig- 
ures 5 and 6) of the dispersive delay line relative to the sweep rate of the function generator, phase 
ripples or fluctuations within the system passband, and amplitude ripples within the system pass- 
band. Each of these departures from the ideal is considered below in terms of classical linear 
system theory. 

Effects of Finite Bandwidth 

The formal derivation of Equation 13 given in Appendix A ignores the fact that the dispersive 
delay line itself, as well as the rest of the transmission system, is a bandpass rather than an all- 
pass structure. The response in the finite bandwidth case can be found using the Fourier integral. 
With reference to Appendix B, the signal x N (t) has a spectrum x N (f ) given by 

x n ( o = e Iff [ Z fel)- Z M ■ f>0 (22) 


where 


and 


a = f Q - f + N5f 


€ 


1 

2 



)• 


The general expression for the finite bandwidth response is then given by 


y N (t) 


/•co 

X N ( f ) H f ( f ) e i2Tr ft df 

J -CO 


(23) 
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where the magnitude of H f (f) is zero except where 


and 


f 0 - T 0 < f < f 0 + 2 mT 0 


f > 0 


H f (-f) = H* ( f ) 

X N<- f ) = X N (O ' 


Then 


y N (O 



rf 0 +Af 

| 

f 0 -Af 




[ z (^l) - z (^ 2 )] ei27rft df 


(24) 


This integral has been evaluated in general terms by DiFranco (Reference 20) and his result 
(for a chirp radar system with r = 50) shown, together with the theoretical | (sin x)/x| curve, in 
Figure 9. Unfortunately, his integrated expression is so complex and lengthy that little is to be 
gained by reproducing it here. Nevertheless, in view of the signal-processing loss and sidelobe 
distortion evident in Figure 9, formulas should be obtained for the peak response and the response 
at the adjacent nulls in accordance with the analytic philosophy outlined above. 

To determine the peak response, we rewrite Equation 24: 


(O 


1 

/2/I 


r f ° +Af r 

[z(0 -Z(£ 2 )l e^'^t+NSf) 

Jf -Af L J 


df 


(25) 



2 M T 0 t 

Figure 9 — Di Franco's finite bandwidth output waveform, T = 50. 
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The peak response is anticipated at t = - NSf/u : 


y N 



l 


r f o +Af 

f 0 -^f 


[z(^)-z(^ 2 )] df • 


(26) 


Introducing a normalized discrete frequency modulation, 8 N = 2N8f/Af , and integrating with 
respect to 6 - a/Af, yields 





Af 




-1+(8N '2) 


f 

) 

(2*-l)J 

M0 


(27) 


Using the formula for the definite integral of Z(x) which is given in Appendix C, we write 
Equation 27 as 


(- = ( ? n + 3) z[/5 (S N * 3)] + (S N - 3) (S N - 3 )]-(S n + l) z[/f (« N + l)] ~ l) z[/J (« N ~ l)]| ' 


'(28) 


This result reduces to the infinite -bandwidth case if F is large enough so that the several 
Fresnel integrals can be replaced by their values at infinity, i.e. 


whence 




r 

8 


I 2 + 2i i 2 


f . 


Cases of interest are $ N = 0, ±1, corresponding to pulse spectra at the center and extreme 
edges of the system pass band. Operation on Equation 28, using asymptotic forms for Fresnel 
integrals of argument greater than 20 (Appendix C), shows that peak response is only down by some 
0.1 db in the cases § N = 0, ±1. The situation will be similar for intermediate values of $ N , and, in 
view of other anticipated system losses, finite bandwidth is not a serious consideration— at least 
as regards losses in processing gain. It remains to be seen how much sidelobe distortion (as is 
evident in DiFranco T s F = 50 case) will be present in the output waveforms. 

The Nth null of the Nth output pulse occurs at t = - 5f/p. (N M) , and we proceed, as above, to 
evaluate Equation 25 for this particular value of time. Again we use the normalized discrete 
frequency modulation S N , as defined above, but this time the integration is with respect to the argu- 
ments of the Fresnel integrals. After undergoing several elementary manipulations, the expression 
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for the output signal amplitude at the Mth anticipated null of the Nth output waveform is 


'n [" 7T ( N + M >] = 


r v^(S N+ 3 ) 

J/_ z{^i) 


, r/ryin N -i) . — 

-i27T(N+M)/2/r^, , \ -i2"(N+M)/2/rf 2 

' 1 _ J . — — Z (^ 2 ) e 

J ' / V*( S N-3) 




( 29 ) 


These integrals may be recognized as truncated Fourier transforms of the Fresnel integral 
from f l9 space respectively into the discrete (N + M) i2/V space. For any combination of the 
indices N and M, Equation 29 can be written out directly in terms of Fresnel integrals using the 
general formula for 


Z(x) e l27Tyx dx , 

given in Appendix C. Because of its length, the integrated expression will not be given explicitly 
here. For the system in question, several computations were made (in particular, for N = 0 and 
M = 1,3; the sidelobe level at these points was more then 50 db below the main response. 

For large I , the limits of the integrals diverge and the integrals tend to cancel. This agrees 
with the result that signal-processing gain is the same for both finite and infinite bandwidth situa- 
tions for large I , and reflects the fact that for large r the linear FM pulse tends to have a rectangu- 
lar frequency spectrum (Appendix B). 

Effects of Slope Mismatch 

With reference to Figures 5 and 6, the dispersive delay line should have a linear group-delay 
characteristic whose slope is the reciprocal of the frequency-vs-time slope of the function generator 
waveform f(t). Departure from exact matching of the overall slopes produces waveform distortion 
in the y N (t). Slope mismatch ("quadratic phase distortion") is discussed in some detail in the 
literature of chirp radar, but it seems that there are no explicit equations for losses in f and the 
behavior of anticipated nulls. The desired equations are derived here by an approximation evolving 
in Reference 19. 

Figure 10 illustrates the desired quadratic phase characteristic of the dispersive delay line, 
and also a phase characteristic corresponding to slope mismatch. With reference to Equation 11 
(here renumbered), 


y N (t) = x N (t)*h(t), ( 30 ) 

the response in the presence of quadratic phase distortion of the type shown in Figure 10 is given 
by the magnitude of the convolution of y N (t) with the impulse response of a fictitious linear filter 
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Figure 10— Quadratic phase distortion. 


of system function 


= e U k -D^-'») 2 . 


(31) 


For slope mismatch the same across the band, the waveshapes of all output signals are the 
same for all N. In the case of N = 0, 


(t) 


f" fy 0 (r)e i t <k ' 1 > ^ f - f » ) 2 U-^)dfdr 

J —00 —CO 


(32) 


The quadratic phase distortion is conveniently characterized by a phase shift <t> measured at 
f 0 ±mT 0 , as shown in Figure 10: 


4> = (k - 1)77 M T 0 2 


(33) 


and, with this convention, Equation 32 becomes 


y* CO 



y 0 ( T ) 




e i27rf(t-r) 


df dr • 


(34) 


It does not appear possible to evaluate this integral exactly. An approximation discussed in 
Reference 19 ignores a residual frequency in y 0 (t) . This is permissible because the effective 
duration of y ( t ) is short. This approximation is equivalent to assuming that the envelope of 
x N (t) is not exactly rectangular but approximately so, like the spectrum of X N (f)(see Appendix B). 
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When this assumption is made, the Fourier transform of y 0 (t) is zero except in rectangular 
bands of width 2 p T 0 centered about ±f 0 . Replacing the time -domain convolution in Equation 32 with 
frequency-domain multiplication and changing variables leads to 


y^(t) - 


ft 


" T ° »<!''* of e i™. df 




(35) 


This expression can be evaluated, using formulas 7.4.38 and 7. 4. 39 from Reference 21, to yield 


y* (t) 


\ 


24> 


^ (rr M T 0 t +<£) 


- Z 


77<f ("MTq t -0) 


Evaluating this result at t = 0 gives a concise formula for the loss in processing gain: 


y*(°> = /w } 


2 d> 

z\V-£ 


(36) 


(37) 


In the limit for small </>, this reduces to the matched case, for 


1 im 


Z(x) 


1 . 


Equation 37 is plotted in Figure 11 as peak loss in db vs 4 , in radians. 

Equation 36, at the time of the Mth expected null, gives the signal amplitude due to slope mis- 
match as 




/MSf\ _ 1 rr—, r 77 [ /~ 2 ~ /Wr \ rr (Wj \ 

nr) = 2 r^- T o r 2$ z \y^ \~r + < *>)J _ z \y 5$ It -*) 


(38) 


Figures 33 and 36 of Reference 19 suggest that special attention be given to the case M = ±1, 
i.e. the leading and lagging anticipated nulls. Figure 12 of this paper plots Equation 38 in normalized 
form (i.e., for an undistorted peak response of unity). 

Equations 37 and 38 and Figures 11 and 12 apply, of course, to all 127 signaling waveforms, 
not merely the N = 0 case. Further, they apply to all rectangular-envelope linear FM pulse- 



Figure 1 1 — Loss in Tvs quadratic phase distortion <f>. 



$ (radians) 


Figure 12— Normalized amplitude of leading 
and lagging nulls vs 4>. 
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compression systems in which the time-bandwidth product r is large enough to justify the approxi- 
mation regarding the residual FM in y N (t ) . 

Effects of Phase Ripple 

Any phase fluctuations of ripples in the transmission system may be considered concentrated 
in the dispersive delay line proper. The nonlinear group delay associated with the system phase 
ripple is given by the general formulation 

, 1 d<p( f ) 

t(f) = - 2^ T * (39) 

An analytic expression for arbitrary phase ripple in the system passband is given by writing 
the delay line system function as 




A(f)e 1«.>.*,<0] , 


where, as before, 


A(f) = 1 

4>(0 = J ( f " f o) 2 

and the phase ripple is expressed as a Fourier sine series* 


(40) 


<£ r (f) 


L 


a. sin j tt 


(tz f o) 
Af 


(41) 


The interval of expansion is 


f 0 - Af < f < f Q + Af , 


where 


Af = 2 u T 0 


is the system semi-bandwidth. 


♦The completness of the representation is not at issue here; the desired phase is neither even nor odd about f Q ; small deviations due to 
tuned circuits, etc., tend to be odd. 
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The response in the presence of phase ripple is given by the convolution of the undistorted 
response y N (t) from Equation 30 with the impulse response of a fictitious filter of system function 


e 


03 

i 2. a . sin 
j=l J 


F(f-f 0 )/Af 


The exponential may be expanded using a well-known Jacobi- Anger formula 


oo 

G ias in/3 = ^ ( a ) e ik/3 _ 

k=“ C0 


(42) 


where j k (a) is the Bessel function of the first kind of order k and argument a. If we discard 
higher-order terms, the impulse response is of the form 


h ph = Jo K) S(t) + }~|j k (a t ) s(t + 2AT ) + f -1 ) k s (t- 2 At)1 


(43) 


where 8(t) is the Dirac impulse function. Convolution with y N (t) produces 


y N (t) * h ph (O = Jo ( a i)y N (O + J2 Jk ( a i) [^n I 1 + 2Af) + ( _ 1 ) k y N ( l _ 2 A 7 )] 


(44) 


where the property 


f(t) * s(t - t 0 ) = f (t - t 0 ) 

of the Dirac function has been used freely. The result in Equation 44 is graphically described as 
’'paired echoes”* and suggests (on the basis of energy conservation) the following addition theorem 
for the squares of Bessel functions: 


j 0 2 (x> + 2 22 h 2 <*) = 1 • 

k=i 

This result was first given by Hansen and is discussed in Section 2.5 of Reference 22. r 

Further discussion of the phase ripple problem is best deferred while amplitude ripple is 
studied. 

*A terminology apparently due to Wheeler (Reference 23)* 

'This may be viewed as a particular consequence of the Parseval theorem. Also, Equation 44 may be considered a time dual of the line 
spectrum of a carrier angle modulated by a sinusoid. 
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Effects of Amplitude Ripple 


Departures from a flat amplitude characteristic within the system passband can be considered 
concentrated in the dispersive delay line proper (assuming that their effects are not altered by 
system nonlinearities such as hard-limiting which may be encountered in satellite repeaters) and 
writing the system function as 


H. (f) = A(f)A a (f)e^< f > , (45) 

where A(f ) and <£(f) have been defined above and A a (f) is expressed as a Fourier cosine* series 

fa (f ~ f 0 ) 

A a( f ) = /_> COS ““ Af • (46) 

fa o 

The response in the presence of amplitude ripple is then calculated by convolving y N (t) with 
the impulse response associated with the system function of Equation 46. Analogous to the phase 
ripple case, the impulse response is of the form 


h a (t) = b o s (t) + \ b-l [ s (t + 2Af) + S ( t - 2&F )_ ' 

■ 1=1 

Convolution withy N (t) again produces paired echoes, but of the form 


(47) 


y N (t) * h a (t) 


b 0 yu( t '> 


i 

+ 2 


b -e [?« (‘ + m) + y N G _ 

i=i 


(48) 


Combined Effects of Phase and Amplitude Ripple 

In view of the basic similarity of the phase-ripple and amplitude-ripple effects, they should be 
combined into a single mathematical formulation. This is accomplished by multiplication of the 
system functions h ph (f) and A a (f) or, equivalently, by convolution of the impulse responses h ph (t) 
and h a (t) . With the latter, the overall impulse response for superposed amplitude and phase 
ripples is of the form 


h ph( t )* b a< t ) ' Jo( a l) b 0 8 ( t ) + ^[jo( a i) T + Jk( a i) b o] S ( t+ ®F) + ^^o( a i) 

k=l k=l 


;+ (-i) k j k K)b 0 




(49) 


♦Again, the completeness of the representation is not at issue. The desired amplitude characteristic is even about the center frequency, 
and small departures tend to be even. 
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If we take only the initial terms (corresponding to the leading and lagging echoes), convolution 
with y N (t) yields 


y N (*)*»■>,,,, oosh.oo = J 0 ( a i) b 0 y N (t)+ J 0 ( a i) ~z + Ji( a i) b o y N ( t+ w) + Jo( a i) t _ Ji KK y N (* ■ sr)-( 5 °) 


Figure 13 illustrates the main response and leading and lagging echoes is given by Equation 50. 
Higher-order terms which have been discarded in the above formulation correspond to echoes 
farther removed from the main response, but these in general will be much smaller than the main 
echoes shown in Figure 13. In the case of arbitrary amplitude and phase ripples, the echo struc- 
ture will not tend toward symmetry; under certain conditions, e.g. 


D i 

Jo ( a i) T = ■ Ji ( a i) b o 


one echo may vanish completely. Further, Equation 50 describes rf signals prior to envelope de- 
tection (and not the envelope functions themselves, as indicated in Figure 13 for simplicity). Hence, 
superposition of leading and lagging echoes on the main response will not necessarily be phase- 
coherent; an rms superposition of the envelope functions would probably be appropriate.* 

Regardless of phase coherence, the loss in r due to paired echoes is approximately 

2 °l°g 10 [jo KKjdb , 

which is plotted in Figure 14 for the case b 0 =1. 



2pT 0 t 

Figure 1 3— Simplified geometry of paired echoes. Figure 14— Loss in T due to paired echoes. 

*Reference 19 contains an exhaustive treatment of the paired-echo superposition problem relative to sidelobe reduction in chirp radar 
systems. It is noted that the phase coherence of leading and lagging echoes was not taken into consideration in that analysis. 
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SUMMARY 


The theoretical operation of a proposed maximum-likelihood Darlington detector for discretely 
coded orthogonal PFM signals of unknown phase has been described. It was shown that in principle 
the communication technique is an ideal M-ary system in the sense defined by Turin. 

Anticipated departures from the theoretical system are identified: finite bandwidth in the 
Fourier transforming dispersive delay line, slope mismatch or quadratic phase distortion, ampli- 
tude ripples, and phase ripples. The departure from the ideal characteristic is represented in 
terms of classical time-invariant linear system theory, and characterized by a system function 
and/or an impulse response. The analysis derives general results for the effects of the several 
departures from the ideal (approximating when necessary). 

Finite bandwidth produces negligible losses in system-processing gain (because of the large 
time-bandwidth product of 128) and little significant sidelobe distortion. Slope mismatch leads to 
serious deteriorations in processing gain and sidelobe waveshape; relatively little quadratic phase 
distortion can be tolerated in implementing a practical system. Amplitude and phase ripples pro- 
duce paired echoes that reduce the processing gain and spread out the main response in time, lead- 
ing to potential intersymbol interference or loss of signaling waveform orthogonality. 

The results presented can be used to set upper limits on the departures from the theoretical 
system, tolerable in implementing the communication technique. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, July 3, 1967 
160-44-05-14-51 
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Appendix A 


The Quadratic Phase Network as a Fourier Transformer; Derivation of Equation 13 


This appendix performs three functions. First, it gives a detailed derivation of Equation 13; 
second, it shows that the pulse-compression phenomenon studied here is only one instance of a 
much more general Fourier transforming property of quadratic phase networks; and third, it pre- 
sents all the required computations in step-by-step form, unlike the main text in which few details 
of the mathematical manipulations were given. 

We wish to compute 


y N ( t ) " y^( t ) ~ x N ( t )* h ( t ) 


where 


X N (t ) - COS 277 


(f 0 +NSf)t 
I t I > T„ 


A*' 


- T < t < T 
A o __ t _ i 0 


and 


h(t) 


£ 


r [( f - f o) 2 /^- ft o] 


i 2-rrf t 


df 


We express x N (t) as the real part of a complex signal x N (t) and solve for y N (t) by taking the 
real part of the complex response. We have 


^(t) = e i2 d( f o +NSf ) t+ ^ 2/2 ] 

= 0 1 1 1 > T 0 . 


-T < t < T 


The spectrum of ^ (t) is given by 


\( f > 



[(f 0 +NSf)r+MT 2 /2] 


e -i277fT 


dr 


25 



and 


/•OJ 

*' -a 


x N (t)*h(t) = (f)H(f)e* 217f 1 df 


ri 


a i 2-n [(f o+N 8 f - f )r+/xr V 2 ] 


dr 


a i2TT[(f-f 0 ) 2 /2M-f(t-t 0 )] 


df 


-r*r-‘ 


( w /^)[(f-f 0 ) 2 + (MT) 2 +2M(f 0 +N8f)r-2MfT+2Mf(t“t 0 )] 


df 


Taking t' = t-t 0 , we rearrange the quantity in brackets to read, after completing the square 
in the first few terms, 


(f + /xt ' - /at- f Q ) 2 + f Q 2 + 2/x(f 0 +NSf)r + (m 7 " ) 2 - (^t' -/xr - f Q ) ; 


This rearranges after several cancellations to read 


(f + /xt'-/xr-f 0 ) 2 + (2/zNSf + 2/x 2 t # )r ~ (/xt ' ) 2 + 2/xf 0 t ' 


We can now write the complex response as separate f -domain and r-domain integrals: 


^ i2ir \f n t'~ 

x N (t)*h(t) - e 


t'-(Mt') 2 ] r ° ei2 , T(NSf+flt . )T ^ r e i 

J -X j - CO 


( 77 /M)( f 0 ) 


df . 


The improper f-domain integral can be evaluated by changing variables according to 


- /tT ( f ~ f o) 


dx 


/T 

f M 


df 


to obtain the tabulated integral 


/•CD 

J e ix 


dx “ J 771 . 
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This yields 


*n( 0 *h(t) 


i 27 r[f t '- 0 tt' 2 / 2 )+l/ 8 ] f ° 

e I 

“ T 0 


e i2^(N8f+Mt # )r ^ 


Integrating over r space gives 


* N (O * h ( t ) 




i 2 w[f # t-(m'*/ 2 )*l/ 8 ] 


sin 2 ttT 0 (NSf + /zt') 
2 tt T 0 (NSf 


The real part of the complex response is 


/ N (t) “ /V T 0 2 COS 277 ^f Q 


M t' 2 i\ Sin 2ttT 0 (NSf + M t^) 


t - 


2 + 8/ 2tt T 0 (NSf +^t T ) 


In regard to the envelope of this expression, disregarding the frequency modulated carrier: 


y N (t) = /VT 0 2 


sin 2 ttT 0 (NSf + /Lit ' ) 
2 ttT 0 (NSf + 


Finally, t' - t-t 0 gives: 


y N (t) 


A /t o 


sin 2 t 7T 0 [/i(t - t 0 ) + NS fj 
2 ttT 0 [ M (t- t 0 ) +NSf] 


which is Equation 13 of the main text. 

The above result is a particular example of a more general phenomenon, namely, the Fourier 
transforming property of the quadratic phase network. Consider an arbitrary signal of real en- 
velope e(t) modulated by a complex chirp or linear FM modulation, i.e. 


V2)] 

x(t) - e(t)e Lo J . 

This chirp modulation may, of course, be superposed on any signal by mixing techniques. Let 
X(f ) be the Fourier transform of x(t), and E(f) be the Fourier transform of e(t ). Now let x(t) 
drive an all-pass quadratic phase network of system function 

H(f) = e^^-'o ) 2 . 
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The spectrum of the response y(t) is 


Y(f) 


H(f)X(f) 


f H( f ) x( t ) e _i217ft dt 


, f €( r ) e i< ^ ) ( f - f o) 2 e -[( f o- f X + W 2 / 2 )] ^ 

* -00 

The integrand may be written with a perfect square in the exponent, 

f i(ir //!)[«- f , 

Y(f) = e(r) e L< - 0) J dr 

» -00 

and the response in the time domain is 


y(t) 


f f e(T)e i< 7 r/M, k f - f »)H 2 e i27Tft df dr . 


Equation 



77 

7^ 


( f ' -MT) 2 


+ yU77t 2 + 277 f' t “ 277/xtT , 


where 


f 7 


f 


f 


o 


gives 


y(t) 


5 i[ 2 rf 0 t-^t 2 ] f" €(T)ei27r/itT f 

j -00 J -00 


df' dr 


TaMng 


we have 


£ = 


( f ' "M-r) 


+ y^77 t 
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and the response is 


y(t) 



€(t) e i2,r ^‘ T d t . 


The remaining r space integral is the Fourier transform of s(r ) into -/it space, 


y(t) = E(-Mt) 

The above result may be summarized as follows: the quadratic phase network is a Fourier 
transformer for the envelopes of chirp- modulated signals. 

In the system studied in the present essay, 



i - T 0 < t < T 0 

e(t) - 

0 I t I > T 0 , 


and its Fourier transform is 


E(f) 



€(t) e' i27,ft 


dt 



sin 277 fT 0 
2T 0 277 ff^~ ' 


Then 


E(-Mt) 



sin 277 T q fJLt 
277T 0/ It~ 


and 


y(t) 


i^T 0 2 


i277[f 0 t-(ftt 2 /2) + l/8] S * n 277T 0 ^it 
e 2tt1 0 ijX 


The envelope of the response is 


y(t) 




sin 277T 0 yu-t 
2 77 T 0 77t 


which is in essence Equation 13 of the main text. 
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Appendix B 


Signal Spectra and System Bandwidth Considerations 

The spectra of the actual signaling waveforms s N (t) are sin x/x centered at ±(f 0 +N 8 f ) and 
are unimportant.* This discussion will be concerned with the spectra of the detector waveforms 
x N (t). It is these spectra that establish the bandwidth requirements for PFM-to-PPM conversion 
or spectrum analysis in the dispersive-delay-line demodulator. Formally, the spectrum of any 
x N ( t ) is given by 


X N (f) = x N (t) e" i27Tft dt 


- I ! (O + I 2 (O 


where 


MO 



^[( f 0 “ f + N 5 f)t + Ci t 2 /2 )] 


rr[(f 0 + f+NSf)t+Git V 2)] 


dt 


dt . 


In each case, the integration is attacked by completing the square in the exponent along the 
lines discussed in Appendix A. After the required manipulations, 



e -i("a 2 /u) 

(f) 

2/2^ 


e -i(^ 2 / M ) 

I 2 <f) 

2/2^ 


[Z(b)-Z(a)] 


[Z(d)-Z(c)] 


*Spectrum truncation due to signal prefiltering is neglected here; as noted in the main text, it is desirable to control the system noise 
bandwidth by such an operation, and this will reduce the bandwidth requirements of the dispersive delay line accordingly. The present 
discussion is based on the use of rectangular-envelope chirp signals for simplicity. 
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where 


(f 0 -f +NSf) 


(f 0 + f +N8f) 



o 

II 

sji 

■£(£*.)= 

. . Jr m 
d - riVsr 


What can be said of these spectra? First, for positive frequencies, X N (f) is practically de- 
termined by I x (f) alone. The arguments c and d are so large for all values of f of interest that 
l 2 (f) is vanishingly small as long as system rf and if bandwidths are small compared to rf and if 
frequencies. This will hold in any practical system. For negative frequencies, i 2 (f) effectively 
defines the spectrum; by symmetry, the signal power spectrum is even about the origin. Hence, 
system bandwidth requirements can be determined satisfactorily from l 1 (f ) . 

Reference 19 shows that the magnitude of i x (f ) is well described by a rectangular function of 
width 2 fj T 0 centered at f 0 + NSf . For large values of the time -bandwidth product P, the approxima- 
tion to a true rectangle is excellent; this is shown in Figures 4, 5, and 6 of Reference 19. 

In regard to the several spectra corresponding to all values of N, virtually all the signal power 
is concentrated in a band of width 4/x T 0 centered about f 0 . This is the required bandwidth of the 
dispersive delay line for effective PFM-to-PPM conversion of unprefiltered rectangular envelope 
system waveforms. 
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Appendix C 


Some Notes on the Complex Fresnel Integral 


The complex Fresnel integral has been used extensively in this study of waveform distortion 
in the Darlington demodulator. This appendix gives definitions, summarizes important properties, 
and presents several useful results pertaining to the integration of Z(x) . 

The complex Fresnel integral is defined by 

Z(x) » I" e l <’"' 2 >^ 2 df . 

Jo 

This holds for all x on the real line and can be extended to complex x as well. Z(x) is a single- 
valued function that possesses derivatives of all orders everywhere. Figure Cl, reproduced from 
Reference 24, shows Z(x) for small arguments in terms of its real and imaginary parts C(x) and 
S(x), i.e. 


Z(x) 

= C(x) + 

iS(x) 



r 



C(x) 

I cos 

(2 H 

df 


Jo 




r 



S(x) 

1 sin 

Jo 

U n 



Also shown in Figure Cl is the fact that if C(x) and S(x) are plotted in rectangular coordinates 
x, C, and s, they appear as the projections of a spiral in space whose projection on the c, S plane 
is the famous spiral of Cornu. 

Important properties of Z(x) are 


Z(-x) = - Z(x) 

Z(0) = 0 

Z(±co) = ± (^±) • 
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At several points, this study has used tables of Z(x) in co input ations and plotting curves. The 
principal source has been Reference 25; asymptotic forms for large arguments were taken from 
Reference 21. 

Some operations on the Fresnel integral using integration by parts are listed below. The 
Leibnitz rule for differentiation of the definite integral is also used: 


_d_ 

dx 



f(f. x)df 


r b(x > d 

f [ b ( x ), x] b ' (x) - f[a(x), x]a'(x) + ^ [f(£. x)] dx . 

J a(x) 


Result 1— Definite Integral of Z(x) 


f Z(x) dx - x 2 Z(x 2 ) - Xj Z(x,) + ^- [e‘ ( /2>X2 - e‘ C /2>X ‘ ] 

•'x . 
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Result 2— Definite Integral of Z(ax + b) 


f Z(ax + b) dx - i [(ax 2 + b) z(ax 2 +b) - (axj +b) z(ax x + b)l + •— [e l(7T/2) ^ aX2+b ^ - e 1< ” /2)< ' aXl+b ’ ) ] 

Jx , 


Result 3— Truncated Fourier Transform of Z(x) 


■L 


F(y ) = Z(x) e~ i27T y* dx 


"2Sy [e— 2 Z(x 2 ) - Z(xj] - [ Z (x 2 - 2y )■ - Z (x, - 2y)] 


Result 4— Fourier Transform of a Function Defined by Z(x) 

f(x) -- z(x +X 0 ) - z(x -x 0 ) 


F(y) 


I f(x) e' 12 "** dx = 2x 0 /2 i e' i2l » 


sin 2rr x n y 
2 ttx 0 y 


Result 3 reduces to Result 2 at the origin, where one must use L T Hospital T s rule to evaluate 
the indeterminate form. Result 4 is of interest, in that it shows that independent of the details of 
the waveshape of a function defined in terms of Fresnel integrals as written above, its Fourier 
transform always has a (sin x)/x envelope. This matter is discussed in connection with passive 
generation of chirp-radar signals in Reference 19. 
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